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Abstract. We present measurements of spin transport in ultracold gases of fermionic 
6 Li in a mixture of two spin states at a Fcshbach resonance. In particular, we study 
the spin dipole mode, where the two spin components are displaced from each other 
against a harmonic restoring force. We prepare a highly-imbalanccd, or polaronic, spin 
mixture with a spin dipole excitation and observe strong, unitarity limited damping 
of the spin dipole mode. In gases with small spin imbalance, below the Pauli limit for 
superfluidity, we observe strongly damped spin flow despite the presence of a superfluid 
core. 
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1. Introduction 

The quality of transport is one of the most important properties distinguishing states of 
matter. Of great technical importance, electrons in condensed matter materials can flow 
as currents or supercurrents, or be localized in an insulator, or even switch their state of 
conductivity through controllable parameters like an applied magnetic field. The task 
of many-body physics is to develop models that may explain the observed transport 
properties in a system. Dilute atomic gases cooled to quantum degeneracy provide ideal 
systems for testing many-body theories. In particular, Feshbach resonances pQ in atomic 
Fermi gases allow experimental control over the strength of two-body interactions, 
giving access to the BEC-BCS crossover regime [2] Transport properties have 
played an important role in characterizing strongly interacting Fermi gases in the BEC- 
BCS crossover, with the observation of hydrodynamic flow indicating nearly perfect 
fluidity jH the measurement of collective excitation frequencies probing the equation 
of state 0, [S] , and the observation of vortex lattices in rotating gases demonstrating 
superfluidity [9]. However, there have been few studies of the transport of spin in 
strongly-interacting Fermi gases. Spin transport has been used in weakly-interacting 
systems to probe the onset of quantum degeneracy [ID] and to generate spin waves [TT] . 

Here we study spin transport in strongly interacting two-component Fermi gases. 
Spin currents are strongly damped in such systems due to the high collision rate between 
opposite spin atoms: as two-body scattering does not conserve relative momentum, each 
scattering event on average reduces the net spin current [T2] . At the Feshbach resonance, 
scattering is maximal, with a mean free path between collisions of opposite spins that 
can be as short as one interparticle spacing - the smallest possible in a three-dimensional 
gas. Recently, we reported measurements of spin transport in strongly interacting Fermi 
gases with an equal number of atoms in two spin states [13]. We demonstrated that 
interactions can be strong enough to reverse spin currents, with two clouds of opposite 
spin almost perfectly repelling each other. The spin diffusivity was found to reach a 
lower limit on the order of h/m at unitarity, the quantum limit of diffusion. Here 
we consider the case where the number of atoms in the two states is unequal, and 
study spin transport in the polaron and phase-separated superfluid regimes. In highly 
polarized systems that remain non-superfluid down to zero temperature [HI [15], [16] , spin 
currents are expected to become undamped due to Pauli blocking [TT] [181 US]- in this 
imbalanced regime, a high-frequency mode observed after a compressional excitation was 
interpreted as a weakly damped spin quadrupole (or breathing) mode [20]. The question 
of the damping properties of the spin excitation and its temperature dependence was 
left open. Spin transport properties of ultracold Fermi gases have been investigated 
theoretically in [21], [22] [TT] US [191 [22] , allowing direct comparison between theory and 
experiment. 

In section [2] we present measurements of the damping rate of spin excitations in 
highly polarized Fermi g function of temperature. We show that damping is 

maximal at finite temperatures. In section [3] we study smaller spin polarizations, below 
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the Pauli limit of superfluidity [Hj, just enough to reveal the presence of a superfluid 
core in the system. We show that the spin dipole mode is strongly damped despite the 
presence of the superfluid. At first this is surprising, as one might expect the superfluid 
to flow without friction, and we discuss an explanation for this observation. 

2. Highly Imbalanced Fermi Gases 

Fermi gases with resonant interactions can remain normal down to zero temperature 
if the spin imbalance exceeds the Pauli (or Clogston-Chandrasekhar) limit [2H 125] HU 
[T5| 120] . We refer to the spin state with the larger population of atoms as the majority, 
or spin up state, and the state with fewer atoms as the minority, or spin down state. 
Radio-frequency (RF) spectroscopy [16] on such systems confirms the quasi-particle 
picture [23 123 EH] where minority atoms are dressed by the majority Fermi sea, forming 
a quasi-particle known as the Fermi polaron. The energy of a single polaron in a 
zero-temperature Fermi sea of spin up atoms has been described using the effective 
Hamiltonian [291 \T7\ [18] 

where p is the momentum of the polaron, m* is the polaron effective mass, fi^ 
is the local spin up chemical potential, and a characterizes the polaron binding 
energy. The parameters a and m*/m, where m is the bare mass of spin up and spin 
down fermions, have been measured experimentally [201 EE I2H1 EI] and calculated 
theoretically [291 E2] E3J EI], giving a = 0.62 and m*/m « 1.2 at zero temperature. 

We consider a mixture of spin up fermions and N± spin down fermions at 
temperature T with equal masses and resonant interactions in a spin-independent 
potential of the form 



where r 2 = x 2 + y 2 . The minority cloud is initially displaced by a small amount Z±(0) 
along the z axis and is allowed to relax to its equilibrium position. 

In the limit iVj. <C iV-j-, the motion of the spin up cloud due to momentum absorbed 
from the spin down cloud may be neglected. The equation of motion of the spin down 
center of mass Z± is then [T7] 



where the factor of (1 + a) is due to the attraction of the minority fermions to the 
majority cloud, and 1/rp is the momentum relaxation rate due to collisions and is 
equivalent to the spin drag coefficient [121 121] • By dimensional analysis, h/r P must 
be given by the majority Fermi energy times a universal dimensionless function of the 
reduced temperature T/T F ^ and the ratio T F ^/T F ^ of the Fermi temperatures, where 
Tp|-(|) is the majority (minority) Fermi temperature. The first two terms in (J3| follow 
from Q, ([2|, and the local density approximation, while the third term is due to 



V(r, z) = -mu 2 z 2 + V r (r), 



(2) 



m*Z i + (1 + a)mu 2 Z i + m*Z i /T P = 0, 



(3) 
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damping and is not captured in (JT1). In Eq. ^ it is assumed that the minority cloud 
is small enough in size compared to the majority cloud that variation in 1/rp with 
density may be neglected. Also, ^ neglects a possible back-action of the minority on 
the majority atoms that might deform the majority density profile. 

In our experimental realization of this transport problem we use a gas of ultracold 
fermionic 6 Li atoms. The 6 Li atoms are cooled sympathetically with 23 Na [35] and loaded 
into a hybrid optical and magnetic trap with an adjustable bias magnetic field [36] . The 
magnetic field curvature provides essentially perfect harmonic confinement along the 
axial (z) direction, while the optical dipole trap (laser wavelength 1064 nm, waist 115 
/iin) provides trapping in the radial directions, with negligible contribution to the axial 
confinement. With this system we perform a collection of time series measurements. In 
each time series we prepare the system in a chosen initial state and observe its evolution. 

At the Feshbach resonance at 834 G, the magnetic moments of "spin up" and "spin 
down" atoms, the two lowest hyperfine states of 6 Li, are equal to 1 part in 1000, as 
their electron spin is in fact aligned with the magnetic field. Inducing a spin current is 
therefore extremely challenging on resonance. However, at lower fields, their magnetic 
moments differ, allowing separation of the two gas clouds by a magnetic field gradient. 
Our experimental procedure for producing these separated clouds is as follows. 

We prepare the system starting with about 1 x 10 7 atoms of 6 Li in the lowest 
hyperfine state, at a total magnetic field of 300 G. A small fraction of atoms are 
transferred to the second lowest hyperfine state using a RF Landau-Zener sweep. The 
mixture is then evaporatively cooled for a variable amount of time by lowering the depth 
of the optical dipole trap from ks x 7/xK to a variable final depth between ks x 0.5/iK 
and ks x 1/iK, where fee is the Boltzmann constant. The optical dipole trap depth is 
then raised to ks x 6/iK, where the zero temperature Fermi energy in the majority state 
is between ks x 0.8/xK and ks x 1.3/iK. 

After the spin mixture is prepared at 300 G, the total magnetic field is reduced 
gradually over 500 ms to 50 G, where the ratio of the magnetic moments of the two 
states is 2.5 and interactions are very weak. A magnetic field gradient is applied along 
the z direction for about 4 ms, imparting a linear momentum of the same sign but a 
different magnitude to each spin state. The clouds are then allowed to evolve for about 
30 ms, and they execute about half of an oscillation period at different amplitudes and 
frequencies (the frequency ratio is 1.6 between spin up and down). When the clouds 
have returned to the center of the trap, their centers of mass are displaced from each 
other by about 200 jum (for comparison, the 1/e radius of the majority cloud in the z 
direction is between 200 and 300 /im at this point). A second gradient pulse is applied 
along the same direction to remove the relative velocity of the two clouds. The second 
pulse also removes most of the total center of mass motion. The total magnetic field is 
then ramped to the Feshbach resonance at 834 G in about 5 ms. At resonance, the two 
spin states have identical trapping frequencies of 22.8 Hz[|J 

| The system as a whole oscillates harmonically along the z direction at 22.8 Hz due to the residual 
center of mass energy. This motion does not affect the dynamics in the total center of mass frame 
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Figure 1. Measuring the spin dipole mode of a highly-polarized Fermi gas. (a) 
and (b) show two-dimensional column density images of the minority and majority 
spin state, respectively, obtained using resonant absorption imaging in one run of the 
experiment. The imaging pulses are each 4 fj,s in duration and separated by 6 fj,s. 
The distance between the centers of mass in (a) and (b) is 34 /im. (c) Density of 
the majority (red circles) and minority (open blue circles) versus potential energy for 
z < Z^, obtained from the images in (a) and (b). The black line shows a fit to the 
majority density at the edge of the cloud where the minority fraction is 5% or less, 
(d) Displacement d of the minority center of mass relative to the majority center of 
mass as a function of time. This time series includes the run displayed in (a-c). The 
curve shows an exponential fit. 



To reach low temperatures, we apply a variable amount of evaporative cooling by 
lowering the depth of the optical dipole trap after reaching 834 G. The time available for 
evaporative cooling is limited to about 0.4 s by the relaxation time of the spin excitation. 
To reach high temperatures, we prepare a hotter cloud at 300 G and heat the system 
further at 834 G by releasing the atoms from the optical dipole trap and recapturing 
them. The depth of the optical dipole trap is then ramped gradually to a final value 
in 80 ms. The final depth is chosen to keep the number of atoms and the temperature 
approximately constant during the subsequent evolution, and corresponds to an effective 
radial trap frequency ranging from 80 Hz for the low temperature data to 250 Hz for 
the high temperature data. After preparing the system at the chosen temperature and 
with a non-zero spin dipole moment, we are left with typically ~ 4 x 10 5 atoms in 
the majority state and N± 4 x 10 4 atoms in the minority state. We then allow the 
system to evolve for a variable wait time t before measuring the densities of the spin 
up and spin down clouds using resonant absorption imaging. Note that we limit the 
population of the majority cloud to ensure that the central optical density is less than 
2, allowing for accurate density measurements. 

Figs. [l]a,b show typical two-dimensional column densities of the two spin states 
after evaporative cooling on resonance. From the column densities we reconstruct 
the three-dimensional densities n a (r,z) of each state a =t,i using the inverse Abel 

because the trapping potential is harmonic in the z direction, and therefore, according to Kohn's 
theorem, the dynamics in the total center of mass frame are equivalent to the dynamics of a system at 
rest [3"7] . 
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transformation. The temperature of the system is determined by fitting the majority 
density as a function of potential energy to the equation of state of a non-interacting 
Fermi gas [351 EE! ( Fi g-[lj=) : ?H,fg = -A~ 3 C 3 /2 (-e^~ v ^), where A = ^2-Kh 2 /mk B T 
is the thermal de Broglie wavelength, B = 1/ksT, the fit parameters are the chemical 
potential // and the temperature T, ( 3 / 2 is the polylogarithm of order 3/2, and 
V e s = V(r, z — Zj-) is the effective potential energy. The fit is restricted to the outer 
edges of the majority cloud, where n^/n^ < x c . We used a cut-off minority fraction of 
x c = 0.05 for all clouds with T < 0.5Tpf. For some of the data with 0.5 < T '/TVt < 1, %c 
was increased to 0.08 to increase the available signal, while for the data with T > 2T F ^, 
x c was increased to 0.15 for the same reason. These increases in x c should not affect 
the accuracy of the thermometry because the system interacts less strongly at high 
T '/Tpf. This is demonstrated by our spin susceptibility measurements for the balanced 
case in [13] that agree with the compressibility above T/Tp ~ 1, showing the absence of 
spin correlations in this temperature regime. Also, measurements on the density profiles 
of imbalanced clouds [391 EDI EI] indicate that the majority density is not affected by 
the presence of minority atoms at high reduced temperatures. For normalization the 
central densities n a (0) of each species are recorded and used to define the central Fermi 
energies Ep a = H 2 kp a ./2m a , with kp a = (67r 2 n (J (0)) 1 / 3 , = m, and = m*, and 
Fermi temperatures T Fa = E Fa /kB- 

Spin transport is measured by observing the time evolution of the center of mass 
separation d(t) = Z±(t) — Z^(t) (Fig. flji), with Z^(t) the center of mass of the majority 
(minority) cloud along the z axis at time t, determined from a two-dimensional gaussian 
fit to the column density. We find that d relaxes exponentially to zero, corresponding 
to an overdamped spin dipole mode, and fit the evolution to an exponential function 
d(t) = doe~ t/r . The fitted spin dipole relaxation times r can be usefully normalized 
using oo z and E F ^ by referring to equation (J3|. When (w z r) 2 ^> 1, the first term in ^ 
may be neglected. In our measurements, r is at least 100 ms, so (w 2 r) 2 > 200. The 
product uj 2 z t is then related to fundamental properties of the gas via 



m* m 

u;t = - '— — . (4) 

(l + a)r P v ' 

The quantity f = hu^r/Ep^ then provides a dimensionless measure of the relaxation 
time in our experiment, and can be compared to theoretical calculations. 

Figure [2] shows the measured values of the dimensionless relaxation time f as a 
function of the reduced temperature T/T F ^. f increases at low temperatures before 
reaching a maximum of 0. 13(3) Ep^ for T/Tp^ = 0.40(6), and decreases at higher 
temperatures. We interpret the behavior of the relaxation time at low temperatures 
as a consequence of Pauli blocking: as the temperature is lowered significantly below 
the majority Fermi temperature, the phase space available for a minority atom to scatter 
goes to zero. The reduction in f at high temperatures is expected: at high temperatures, 
1/rp is essentially given by the collision rate in the gas [40] . 1/rp ~ no~v. The scattering 
cross section a on resonance for T 3> T F ^ is given by the square of the de Broglie 
wavelength and is thus proportional to 1/T, while the average speed v of the particles 
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Figure 2. Normalized relaxation time of the spin dipole mode of a highly-polarized 
Fermi g function of the reduced temperature T/Tp^. Tp^ is the local Fermi 

temperature at the center of the majority cloud. The solid curve is the low temperature 

limit from p~7], given by equation ( 
The inset shows the average ratio o 



it. 



Ft 



The dashed curve is the expression 0.08a, t 
"the minority cloud size to the majority clouds size 
as a function of the reduced temperature T/Tp^. The cloud sizes are defined as the 
1 /e radii along the z axis, estimated by fitting a two-dimensional gaussian function to 
the column densities of the two spin states. In both figures, each point is a weighted 
average of the results from 1 to 3 time series, with each time series containing on 
average 30 spin up-spin down image pairs. The error bars give standard deviations 
due to statistical fluctuations within a time series. Where the results of more than one 
time series are averaged, the error bars show the standard deviation of the weighted 
mean, determined from the standard deviations from each time series. 



is proportional to a/T. Hence f is expected to decrease like hnav/Ep^ oc ^jTp^/T. 
We observed behavior similar to Fig. [2] in three-dimensional Fermi gases with resonant 
interactions and equal spin populations in [13], although we see more significant Pauli 
blocking here than in [13] at comparable temperatures. 

It would be interesting to have available a calculation of spin transport coefficients 
such as 1/rp at arbitrary temperatures for comparison with our data. A full solution is 
available for Fermi gases with equal populations in one spatial dimension fH\ 122] and 
shows qualitative features similar to our data, showing a maximum of the spin drag 
coefficient (analogous to f) at finite temperatures on the order of T F . 

We expect our data to differ quantitatively from predictions for a homogeneous 
system. The measured quantity f is a global property of the trapped system, while the 
momentum relaxation rate 1/rp is a local quantity. For T ^> Tpj-, 1/rp oc n^, and 1/rp 
increases with increasing majority density, while for T <C Tpf, due to Pauli blocking 
1/rp oc E F ^(T/Tp^) 2 oc n^ 2 ^ 3 , and 1/rp decreases with increasing majority density. 
Additionally, the variation of 1/rp should cause the spin current to be non-uniform. 
The effect of inhomogeneity should be greater at high reduced temperatures, where the 
minority cloud size approaches the majority cloud size. The inset in Fig. [2j shows the 
ratio of the cloud sizes R^/R^ as a function of the reduced temperature, where R\(i) 
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is the 1/e width in the z direction from a two-dimensional gaussian fit to the majority 
(minority) column density. Indeed R^/R^ increases with increasing T/Tp^. Even at 
the lowest temperatures, R^/R^ remains significant, attaining a value of 0.7, due to the 
finite minority fraction N^/N^ ~ 0.1. The effect of inhomogeneity is therefore reduced 
at low temperatures, but should remain present. 

We compare our results at low temperatures to the low temperature limit in 
Ref. [T7], which can be written as 



a I m \ I 1 



T 



1 + a \ m J \Tj 



■Ft 



(5) 



for temperatures T <C Tp|- [§| The prefactor c changes slightly from c = = 6.89 . . . 
to c ~ 6.0 as the temperature rises from far below Tp±, where even the minority cloud 
is degenerate, to temperatures where Tp\, «C T <C Tp^ and the minority is a classical 
gas [IT]. In our coldest data, T w 0.5TW and Tp± ~ 0.37Vj-, assuming m* = 1.2m. 
To compare our data to Ref. [TT] using (p) we set c = a = 0.6 and m* = 1.2m. 
The experimental data agree with the theory at the lowest temperatures measured 
(see Fig. [2|. The deviation at higher temperatures is expected as the T <C T F ^ limit 
becomes inapplicable. The convergence of the experimental data to the theoretical 
value at low temperature despite the inhomogeneity of the system may be partly due 
to the reduced minority cloud size at low temperatures, which reduces the effects 
of inhomogeneity, as discussed above. Additionally, the variation of the momentum 
relaxation rate with density will to some extent cancel at moderately low temperatures, 
as 1/rp changes from increasing with increasing density at high reduced temperatures to 
decreasing with increasing density due to Pauli blocking at low reduced temperatures. 
The crossing of the experimental curve with the predictions for a uniform system at low 
temperatures therefore does not necessarily indicate that the inhomogeneity is negligible 
at low temperatures in this measurement. 

At high temperatures T ^> Tp^, the spin transport properties of a trapped system 
can be calculated from the Boltzmann transport equation. For vanishing minority 
fraction, we find (now with a = and m* = m and assuming harmonic trapping in 
all three directions) [13] 



T Ft 0.16 /Tp t 

where e = 1 when the drift velocity distribution is uniform. This result features the 

expected dependence cx \J-ip- on temperature. The relative velocity between the two 
spin states cannot be truly constant in space but has to be depressed in the center, 
where the density is highest and momentum relaxation is fastest. In general, 

/d 3 ™ drift (^)e-^ 

§ We omit a term due to the relative velocity of the spin up and spin down clouds, which produces a 
correction of less than 1% in the overdamped, finite temperature regime accessed in this experiment. 
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where V is the trapping potential (here assumed to be quadratic), and v drift (r) is the 
average velocity along z of minority atoms at position r. For example, for a quadratic 
drift velocity profile, v drift (r) = ax 2 + by 2 + cz 2 , the predicted f is reduced by factor of 
e = 2. We find that the high temperature result (|6p with e = 2 leads to close agreement 
with our experimental results (Fig. [2]). This model is interesting because it estimates 
the effects of inhomogeneous density and velocity distributions, but it has shortcomings. 
The drift velocity should remain non-zero everywhere, rather than going to zero at the 
origin as in the quadratic case, and should have a radial component. A full quantitative 
description of the overdamped spin dipole motion in the high temperature limit in an 
external trapping potential will therefore be more complex. 

3. A Superfluid with Small Spin Polarization 

We extend the method of the previous section to study spin transport in Fermi gases with 
resonant interactions and small spin imbalance. When the global polarization jy-^ is 
less than about 75% in a harmonically trapped Fermi gas at low temperature and with 
resonant interactions, the system phase separates into a superfluid core surrounded by 
a polarized normal state region [TU [TSJ [20]. The superfluid core is visible as a sharp 
reduction in the density difference of the two spin states [15]. The transition between the 
superfluid and the imbalanced normal regions forms a sharp interface below a tricritical 
point, where the density imbalance jumps between the two regions (39]. Scattering and 
spin transport at the interface between a normal and superfluid Fermi gas have been 
considered theoretically in Refs. [HI [TO] . 

To observe spin transport in an imbalanced gas containing a superfluid, we prepare 
a spin mixture with a global polarization of 17(3)%. The gas is cooled at 300 G and again 
at 834 G after creating the spin dipole excitation as described in the previous section. 
Two off-resonant phase contrast images are taken to measure the densities of each spin 
state. An imaging pulse tuned halfway between the resonance frequencies of the two 
states directly measures the difference in the column densities (Fig. |3^i) while a second 
pulse, red-detuned from both states (Fig. [3jo) , provides additional information needed 
to reconstruct the total column density in each state [39]. From the column densities 
of each state we obtain three-dimensional density distributions using the inverse Abel 
transformation. 

The two-dimensional spin density (Fig. |3^i) and three-dimensional spin density 
(Fig. [3^) show a reduction near the center of the trap, with the three-dimensional 
density going to zero, characteristic of the superfluid core in imbalanced Fermi gases [15]. 
We have checked that the shell structure remains even after the spin density reaches 
equilibrium. Additionally, we estimate the temperature T of the system to confirm 
that it is cold enough to contain a superfluid. In unpolarized systems, the superfluid 
transition is predicted to occur at about T c /Tp = 0.173(6) |42j, where fc^Tp = Ep = 
h 2 (6n 2 n(0)) 2 / 3 /2m and n is the density per spin state. This theoretical value agrees 
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Figure 3. Spin dipole mode of an unbalanced Fermi gas with a superfluid core. 
Phase contrast images are taken with imaging light detuned (a) halfway between the 
resonance frequencies of the two states and (b) at large red detuning from both states. 
The image in (a) is proportional to the difference in column densities of the two 
states. The depletion of the density difference in the center of the cloud indicates the 
superfluid region, ft is displaced from the center of the majority due to the spin dipole 
excitation, (c) shows the difference in reconstructed three-dimensional densities of the 
spin up and spin down clouds as a function of the z coordinate for z > 0. The depletion 
in the center again indicates pairing and superfluidity [15] . An elliptical average over a 
narrow range of the radial coordinate r is used to increase the signal to noise ratio, (d) 
Three-dimensional densities of the two states as a function of the effective potential 
defined in the text. Solid red circles: majority density, open blue circles: minority 
density. The curves are fits to the densities using the equation of state of a unitary 
Fermi gas at zero imbalance to get upper and lower bounds on the temperature. 



well with the determination of T C /T F by our group []}] Fitting the equation of state 
of a unitary Fermi gas at zero imbalance |43j to the majority (minority) density gives 
an estimate T^m of the temperature. The fits are restricted to V^ e fr > 0.3 [iK, where 
V^eff = V(r, z — Z a ), to exclude the putative superfluid region. Compared to a balanced 
gas at unitarity with NL = N[ = N-f, and at the same temperature T, the majority 
cloud should have a larger size because the interaction energy between the spin up and 
spin down atoms is attractive. We therefore expect that is an overestimate of T. 
Likewise, we expect Tj. < T, and we consider and Tj. to provide approximate upper 
and lower bounds on T. 

Figure |4^i shows the temperature bounds during the approach to equilibrium. Time- 
averaging gives 0.12(2) < T/T F < 0.15(2), where T F = T F ^ T f ±. The error estimates 
refer to standard deviations. These bounds confirm that the system is in the vicinity 
of the superfluid transition. In contrast to the dissipationless flow that defines the 
superfluid state, however, we observe strong damping of the spin dipole mode. Figure 
[4)3 shows that the displacement d between the majority and minority centers of mass 
along the z axis relaxes gradually to zero, rather than oscillating as would be expected in 
a dissipationless system. The 1/e relaxation time r = 360 ms corresponds to a spin drag 
coefficient |12l [21] of oj 2 z t = 0. 06(1) Ep^/h, close to the maximum spin drag coefficient 



|| Mark Ku, Ariel Sommer, Andre Schirotzek, and Martin W. Zwierlein, in preparation. 



Spin Transport in Polaronic and Superfluid Fermi Gases 



11 




Time (s) Time (s) 

Figure 4. (a) Reduced temperature as a function of time during relaxation of the 
spin dipole excitation in a spin polarized Fermi gas containing a superfluid region. Red 
solid (blue open) circles: T-^^/Tp^^ from a fit to the edge of the majority (minority) 
spin state using the equation of state of an unpolarizcd unitary Fermi gas, giving an 
upper (lower) limit to true temperature, (b) The displacement of the spin up and 
spin down centers of mass relaxes exponentially, indicating strong spin drag despite 
the presence of a superfluid. 



in non-polarized trapped Fermi gases at unitarity [T3] . 

The strong damping of the spin dipole mode, or spin drag, is surprising when one 
considers that the superfluid moves through a background of spin up atoms at a velocity 
~ d/r = 30 /irn/s, much less than the critical velocity of about 0.3vp (SI US], where 
Vp = \Z2Epjm = 50 mm/s. However, unlike in a superconductor, where the electrons 
move through a stationary background of scatterers, here the background is a gas of 
atoms with a root-mean-square ve.ocity of at .east significantly exceeding the 

critical velocity. In the inhomogeneous case realized here there is an additional source of 
damping due to the superfluid-to-normal interface. One then expects that a fraction of 
atoms in the normal region are reflected at the interface due to the pairing gap [HI [T_H] . 
The relative importance of the two effects depends on whether a spin current flows 
through the superfluid or around it. However, we are not able to determine the spatial 
distribution of the spin current with our present data. 

4. Conclusions 

In this work we presented our measurements on the damping of the spin dipole mode in a 
highly polarized Fermi gas with resonant interactions, over a wide range of temperatures. 
The damping is seen to become weaker at temperatures significantly less than the 
majority Fermi energy, as expected from Pauli blocking, i.e. the fact that quasi- 
particles in a Fermi liquid become long lived at sufficiently low temperatures. These 
measurements provide the first quantitative test of theoretical calculations of the spin 
transport properties of highly polarized Fermi gases. We also observe spin transport in 
a Fermi gas with low spin polarization containing a superfluid region. It is found that 
the spin dipole motion remains strongly damped, revealing the importance of friction 
between the superfluid and the normal component, possibly accompanied by reflection 
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processes at the interface. 
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